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Abstract. In this paper the current status of y-ray observations of starburst galaxies from hundreds 
of MeV up to TeV energies with space-based instruments and ground-based Imaging Atmospheric 
Cherenkov Telescopes (lACTs) is summarised. The properties of the high-energy (HE; 1 00 MeV < 
E < lOOGeV) and very-high-energy (VHE; E > lOOGeV) emission of the archetypical starburst 
galaxies M 82 and NGC 253 are discussed and put into context with the HE y-ray emission detected 
from other galaxies that show enhanced star-formation activity such as NGC 4945 and NGC 1068. 
Finally, prospects to study the star-formation - y-ray emission connection from Galactic systems to 
entire galaxies with the forthcoming Cherenkov Telescope Array (CTA) are outlined. 
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INTRODUCTION 

Starburst galaxies have exceptionally high star-formation rates (SFRs) in very localised 

regions (also called starburst region) that have typical sizes of hundreds of parsecs. In 
these regions, the high SFR leads to an increased formation of massive stars and hence 
supernova explosion rate compared to the rest of the galaxy. The high thermal pressure 
in the central regions often leads to the formation of a starburst wind. More than 20 years 
ago, it has been proposed that starburst galaxies could be sources of HE and VHE y-ray 
emission [e.g. 1, 2]. The basic idea is that particles are accelerated in the numerous 
supernova remnant (SNR) shock waves, escape into the starburst region and interact 
with radiation fields, magnetic fields and with the dense gas. Electrons (and positrons) 
predominantly lose energy via inverse Compton and synchrotron processes, whereas 
accelerated protons (and heavier nuclei) predominantly lose energy in proton-proton 
interactions in which they produce neutral and charged pions which can decay into e.g. 
electrons (positrons) and /-rays, respectively. 

The detection of HE and VHE y-ray emission from the archetypical starburst galaxies 
M82 and NGC 253 with the Fermi-LAT, VERITAS and H.E.S.S. instruments [3, 4, 5] 
has triggered a lot of theoretical work. These studies are aiming to explain the origin 
and properties of the y-ray emission in different scenarios. Most calculations [see e.g. 
6, 7, 8, 9] consider diffuse y-ray emission from inelastic proton-proton interactions 
and subsequent ;r'^-decay as leading loss mechanism, although it has recently been 
shown that individual sources such as pulsar wind nebulae (PWNe) could significantly 
contribute to the y-ray signal [10, 11]. 

The first part of this work focuses on the connection between massive star forma- 
tion and y-ray emission in individual Galactic y-ray sources and links these to galaxy- 
size systems such as starburst galaxies. The spatial and spectral properties of the y-ray 



emission from NGC 253 and M 82 and implications for the origin of the emission are 
discussed in more detail in the second part of this paper. In the third section the y-ray 
emission from the starburst galaxies with Seyfert 2 nuclei NGC 4945 and NGC 1068 is 
addressed. Finally, the y-ray emission from other Local Group galaxies is discussed and 
prospects for the Cherenkov Telescope Array (CTA) to investigate the connection be- 
tween massive star formation and high-energy particles on all spatial scales is outlined. 

STAR-FORMATION AND /-RAY EMISSION ON DIFFERENT 

SPATIAL SCALES 

At TeV energies, the population of Galactic VHE y-ray sources clusters tightly along the 
Galactic plane. It has a scale height similar to that of molecular gas and traces regions of 
massive star formation. Although ~ 1 /3 of the source population is still unidentified 
and lacks counterparts at lower energies, the majority of TeV sources is coincident 
with PWNe, SNRs, SNRs interacting with molecular clouds or stellar clusters [e.g. 
12, 13, 14, 15]. At GeV energies, more and more of the objects that show bright TeV 
emission also emerge as HE y-ray sources in the Fermi-LAJ energy band. Especially 
sources coincident with regions of active star formation show y-ray spectra that smoothly 
connect from hundreds of MeV to multi-TeV energies. The most prominent sources, 
where such spectra are observed are SNRs at different evolutionary stages that interact 
with molecular clouds, e.g. W28 [16, 17], W49 [18, 19], W51C [20, 14] or IC443 
[21, 22]. The preferred interpretation of the y-ray emission in these systems is that 
protons interact with the dense gas and produce ;r°-decay y-ray emission. 

The progenitors of massive stars that undergo supernova explosions at the end of their 
lives typically form and evolve in associations or stellar clusters. In stellar clusters, the 
fast winds of massive stars can merge and/or interact with the SNR shock fronts forming 
a collective bubble or even a superbubble. Existing shocks get further amplified and 
particles can get accelerated to very high energies [e.g. 23, 24]. H.E.S.S. observations 
of the massive stellar cluster Westerlund 1 revealed an emission region ~ 20 times the 
size of the cluster [15]. A very modest amount of energy injected by the numerous 
supernovae that exploded over the past millions of years and the collective stellar 
winds can easily explain the emission level. The morphology and spectrum of the 
VHE emission suggests that a significant part of the non-thermal particles originate 
from Westerlund 1, and that protons are likely responsible for the emission. At GeV 
energies, three of the four HE y-ray sources that are coincident with the TeV emission 
are marked as possibly confused with Galactic diffuse emission, indicating that they 
might physically be connected and possibly associated to the H.E.S.S. emission [25]. 
Future studies will show whether or not the GeV and TeV emission is indeed related and 
how particle acceleration, particle escape and interaction with the surrounding material 
of Westerlund 1 works in detail. The recent detection of HE y-ray emission from the 
vicinity of the Cygnus superbubble with the Fermi-L AT between 1 GeV and 100 GeV 
further supports the idea of particle acceleration in stellar clusters and superbubbles [26]. 
Also in this case the morphology and spectrum (that is consistent with an extrapolation 
to the measurement by Milagro [27]) suggests that protons are responsible for the non- 
thermal emission. 



In summary, the study of y-ray emission from SNRs interacting with molecular 
clouds, stellar clusters and superbubbles with Cherenkov telescopes and the Fermi-LAI 
provided strong evidence for proton acceleration in regions of massive star formation. 
All these Galactic systems show y-ray emission on spatial scales from tens of parsecs 
up to ~ 100 parsec. The y-ray emission from starburst nuclei can be used to investigate 
the connection between massive star formation and non-thermal particles on hundreds 
of parsec scales and to study populations of particle accelerators. It is complementary to 
observations of individual Galactic objects and a powerful tool to probe the paradigm of 
cosmic-ray acceleration in SNR shells. 

7-RAY OBSERVATIONS OF STARBURST GALAXIES 

The detection of VHE y-ray emission from the starburst galaxies NGC 253 and M 82 
using ground-based Cherenkov telescopes has been reported in 2009 by the H.E.S.S. [5] 
and VERITAS collaborations [4], respectively. Both galaxies are right at the sensitivity 
limit of current generation lACTs and it required more than 100 hours of good-quality 
data to establish these objects as TeV y-ray emitters. M 82 and NGC 253 are the first 
two external galaxies that have been detected in VHE y-rays where the emission is not 
powered by a supermassive black hole (as opposed to e.g. Active Galactic Nuclei or 
radio galaxies). Shortly after the detection at TeV energies, the Fermz-LAT collaboration 
announced the discovery of these two starburst galaxies at GeV energies [3]. Also the HE 
y-ray signals are very weak and ?a 1 1 months of LAT data were required to detect M 82 
and NGC 253 between a few hundred MeV and a couple of GeV energies. No indication 
of variability at any y-ray energy could be detected from both starbursts, supporting the 
idea that the emission is not triggered by a supermassive black hole at the centres of 
the galaxies. In the following the observational properties of the y-ray emission from 
NGC 253 and M 82 are presented and the implications for the origin of the emission are 
outlined. A more detailed discussion of theoretical implications of the y-ray emission 
from starburst galaxies can be found in these proceedings [29]. 



NGC 253 

Based on 30 months of LAT data and fa 180 hours of good-quality H.E.S.S. data, 
[9] presented a detailed analysis of the y-ray emission from NGC 253. The angular 
extension of NGC 253 of ^ 30' x 7' makes it possible to disentangle the TeV signal 
from the disc of the galaxy and the starburst nucleus. Indeed, the upper limit on the 
extension of the TeV emission of < 2.4'{3o) implies that the VHE y-ray emission is 
consistent with originating in the starburst nucleus. This is also visible in Figure 1, where 
no significant disc component in NGC 253 is seen at TeV energies. Due to the limited 
angular resolution of the LAT, such a discrimination is not possible at GeV energies. 
The best fit position of Fermi is however consistent with the best fit position of H.E.S.S. 
and the optical centre of the galaxy [9]. 

The large data sets gathered in HE and VHE y-rays allowed to study the spectrum of 
the y-ray emission in great detail. The Fermi spectrum between 0.2 GeV and 200 GeV 
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FIGURE 1. H.E.S.S. sky map of NGC 253 (right, [9]). The zoom into the central part shown at the top 
left is a composite image with optical in red (ESO/DSS), near-infrared in green (2MASS, Ks band) and 
VHE 7-rays in blue (H.E.S.S.)- The zoom on the bottom left shows the central molecular zone of NGC 253 
in molecular Une emission [28] (Reproduced by permission of the AAS). 



is best described by a power law in energy with index The = 2.24 ± 0. 14stat ± O.OSsys. 
The H.E.S.S. best fit spectral index above the energy threshold of 220 GeV is Fyhe = 
2.14 ± 0.1 8stat ± O.SOsys and hence compatible within statistical errors with the LAT 
index. Interestingly, a combined fit of the Fermi and H.E.S.S. spectrum over four decades 
in energy results in a photon index of F = 2.34 ± 0.03 and a fit probability of 30%. In 
case the emission is of hadronic origin, and assuming canonical values for the supernova 
explosion energy EsN = 10^^ erg and conversion efficiency to cosmic rays ^SN Eqr = 
10%, the measured y-ray flux suggests that ^ 20 — 30% of the protons interact with the 
dense material in the starburst region before leaving the system in the starburst wind. 
Although the available statistics is quite limited, the smooth alignment of the GeV and 
TeV spectrum is suggestive of the fact that one energy loss mechanism dominates from 
hundreds of MeV up to several TeV. In this interpretation energy-dependent diffusion 
would play a minor role and (energy-independent) escape of particles from the starburst 
nucleus and proton-proton collisions would be the dominant energy loss mechanisms of 
the cosmic-ray population in the starburst of NGC 253. 



M82 

M82 is with an apparent size of ~ 11' x 4' significantly smaller than NGC 253. It 
is hence challenging to discriminate the starburst contribution of the HE and VHE y- 
ray emission detected by the LAT and VERITAS from the disc of the galaxy [3, 4]. 
However, as for NGC 253, the GeV and TeV emission is consistent with originating from 
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FIGURE 2. 7-ray luminosity spectra of SNRs interacting with molecular clouds and the Cygnus 
superbubble (left) and of star-forming galaxies (right). GeV data has been obtained with Fermi, distances 
used to calculate the luminosities are given in [17, 20, 22, 26, 32]. TeV data for W28 is from H.E.S.S. 
[16], for W51 from MAGIC [14], for Cygnus from Milagro [27] and for IC 443 from VERITAS [21]. 



a point-like source located in the starburst core of IVI 82. Based on three years of LAT 
data, the y-ray spectrum of IVI 82 between lOOlVIeV and 100 GeV can be described by 
a power law with spectral index The = 2.2 ± 0. Istat- The measured VERITAS spectrum 
obtained with a total of ~ 140 hours of good quality data has an index of Fyhe = 
2.5 ± 0.6stat ± 0.2sys, and is within large statistical errors consistent with the Fermi index 
and extrapolated Fermi flux. Within errors, also the flux from M 82 at GeV and TeV 
energies is comparable to the flux measured from NGC 253. 



NGC4945 and NGC 1068 

Similar to the two starburst galaxies discussed in the previous chapter, also NGC 4945 
and NGC 1068 are galaxies with circumnuclear starbursts. These two objects however 
also have radio-quiet obscured Active Galactic Nuclei (AGN). In 2010, both galaxies 
were reported to emit HE y-rays [30, 31]. Using three years of LAT data, [32] analysed 
the spectrum in more detail and reconstructed a spectral index of F = 2.1 ±0.2 and 
r = 2.2 ±0.2 for NGC 4945 and NGC 1068, respectively. The spectra as shown in 
Figure 2 (right) are similar to those of IVI 82 and NGC 253 and other Local Group 
galaxies such as M 31 or the Large IVIagellanic Cloud (LIVIC). So far, no TeV emission 
from these two Seyfert 2 galaxies has been reported, although a detection is within reach 
of current generation lACTs. The HE y-ray emission from NGC 4945 and NGC 1068 
could have their origin in the starburst and/or a possible AGN activity. One way to 
disentangle these components is to study the y-ray lightcurve, as the starburst emission 
is expected to be non-variable. Although statistics is limited, the GeV lightcurves of 
NGC 4945 and NGC 1068 show no indication for variability, which could point towards 
a starburst origin of the y-ray emission. A second possibility is to compare the star- 
formation rate and available target material in these systems with the observed y-ray 
luminosity. In the case of NGC 4945, the level of observed y-ray emission is consistent 
with the estimated supernova rate. Together with the non-variable signal this implies 



that the starburst activity alone can account for the observed GeV signal. Although the 
supernova rate in NGC 1068 is comparable to the ones estimated for M 82 and NGC 253 
the observed radio and /-ray fluxes are an order of magnitude higher. This suggests 
that the AGN component might dominate over a possible starburst activity at these 
wavelengths. However, apart from NGC 4945 and NGC 1068, no other Seyfert galaxy 
with radio-quiet AGN was observed to emit HE y-rays, with upper limits up to one 
order of magnitude below the level seen for NGC 1068 [e.g. 33]. Whether the observed 
emission from these two Seyfert galaxies is indeed related to the star- formation process 
or if it is rather driven by AGN activity requires further monitoring by the LAT and 
search for possible variability as well as TeV observations and more detailed modelling 
of the spectral energy distributions. 

7-RAY OBSERVATIONS OF LOCAL GROUP GALAXIES 

Starburst galaxies are not the only galaxies that have been observed to emit y-rays. 
Firstly, the Milky Way itself is a very strong source of diffuse HE y-ray emission, origi- 
nating mainly from cosmic -rays that illuminate the Galactic interstellar medium and lose 
energy via Bremsstrahlung and inverse Compton processes as well as pion-production 
and decay. Also the Milky Way satellites, the LMC [34] and the Small Magellanic Cloud 
(SMC, [35]) have been detected with the LAT. Moreover, the Andromeda galaxy (M 31, 
[36]) has been reported to emit GeV y-rays. The same authors found a simple scaling 
relation of SFR and y-ray luminosity for these Local Group galaxies that also extents 
to NGC 253 and M 82. [32] used a larger 3-year Fermi data set and examined the y-ray 
properties of 69 dwarf, spiral, luminous and ultra-luminous galaxies in GeV y-rays. A 
quasi-linear scaling relation between y-ray luminosity and the SFR as traced by infrared 
and radio continuum emission further supports the relation between massive star forma- 
tion and high-energy particles. Interestingly, there is a spread of one order of magnitude 
in this correlation that could point towards differences in e.g. the efficiency with which 
massive stars via e.g. supemovae transfer energy into non-thermal particles, how high- 
energy particles interact with the surrounding material and/or how they escape without 
interaction. [32] conclude that between 4% and 23% of the intensity of the isotropic 
diffuse component as measured with Fermi could be accounted for by unresolved star- 
forming galaxies between redshift < z < 2.5. Furthermore, up to ten galaxies might be 
detectable within ten years of operation. 

PROSPECTS FOR CTA 

Figure 2 (left) shows the y-ray spectra of three prominent SNRs that are interacting 
with molecular clouds and the spectrum of the superbubble in the Cygnus region. For 
all these sources, hadronic emission scenarios are favoured over leptonic scenarios. 
The same is true for the y-ray spectra and the emission from star-forming galaxies as 
shown in Figure 2 (right). Especially the smooth connection of GeV and TeV spectra 
of individual Galactic objects and the diffuse emission from e.g. NGC 253 and M 82 is 
striking. So far, the correlation between star-formation and y-ray emission has only be 



studied for GeV-detected galaxies. At TeV energies, only the two archetypical starburst 
galaxies have been detected and no emission from M 31, the LMC, SMC or the Milky 
Way has been reported. This situation will change once CTA is in operation [38]. CTA 
will deliver an order of magnitude better sensitivity, a factor ~ 5 better point spread 
function and broader energy coverage and will allow to study individual objects and 
populations of sources in unprecedented detail. First of all, it will be possible to probe 
whether or not the NGC 253 TeV emission is point-like or if it has an extension, similar 
to the spatial extent of the central molecular zone in the starburst region. Secondly, it 
might be possible to detect the disc of NGC 253 in TeV y-rays in a deep exposure, as 
the expected flux is about one order of magnitude lower than the starburst emission 
[e.g. 9]. Furthermore, the y-ray spectrum will be measured with very good accuracy 
and to lower energies, allowing to investigate if spectral features arise that indicate 
a change in the dominant energy-loss mechanism, or if significant signatures for 7-7 
pair production or a dominance of PWNe at TeV energies become apparent. As for 
NGC 253, also the M 82 7-ray spectrum will be studied in great detail and different 
emission scenarios can be tested. All other Local Group galaxies that have been detected 
with Fermi can also be studied with CTA and the SFR - 7-ray luminosity relation can be 
probed with these systems at TeV energies. Finally, it will be possible to study the whole 
young Galactic SNR population with CTA [39]. With the detailed understanding of how 
particle acceleration, escape, and interaction with the surrounding medium works, it 
might be possible to get insights on how the star-formation process is linked to high- 
energy particles. The study of superbubbles with CTA might provide the link between 
individual 7-ray emitting SNRs and the emission from starburst galaxies. 
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